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Abstract 

We perform a NLO QCD analysis of deep-inelastic scattering data, in which we account 
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in 
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(^ , for absorptive corrections. These corrections are determined from a simultaneous analysis 

of diffractive deep-inelastic data. The absorptive effects are found to enhance the size of 

the gluon distribution at small x, such that a negative input gluon distribution at Q^ = 1 

Q^' GeV^ is no longer required. We discuss the problem that the gluon distribution is valence- 

o ■ like at low scales, whereas the sea quark distribution grows with decreasing x. Our study 

^ ■ hints at the possible importance of power corrections for Q^ ~ 1-2 GeV^. 

In a previous paper [1 , we performed a QCD analysis of the high precision data for diffractive 
/S ■ deep-inelastic scattering, recently obtained by the ZEUS El and HI [4| Collaborations at 
d ' HERA. This process may be denoted 7*p -^ Xp, where the (slightly deflected) proton and the 
cluster X of outgoing hadrons are well-separated in rapidity. The large rapidity gap is believed 
to be associated with Pomeron exchange. The diffractive events make up an appreciable fraction 
of all (inclusive) deep-inelastic events, 'j*p — > X. We will refer to the diffractive and inclusive 
processes as DDIS and DIS respectively. A long-standing question concerns the treatment of 
diffractive events in a global parton analysis of DIS and related hard-scattering data. Are 
they already included in the input distributions or must some account be taken of them in 
the DGLAP evolution? The present paper addresses this question. We show that DDIS is 
partially included in the starting distributions and partially must be allowed for in the DGLAP 
evolution. 

The crucial observation is that, in perturbative QCD, the Pomeron singularity is not an 
isolated pole, but a branch cut, in the complex angular momentum plane 0. That is, the 
Pomeron wave function consists of a continuous number of components. Each component i 
has its own size, l//ii. This property was the basis of the analysis of the DDIS data made in 
[Tj. The advantage of describing the DDIS data using an approach where the dependence on 
the scale fi is explicit is the possibility to use the results to evaluate the absorptive corrections 
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Figure 1: The two-Pomeron exchange contribution to F2. The equahty shows the apphcation of 
the AGK cutting rules, and the relative magnitudes of the cut diagrams. All the permutations 
of the two gluon ladders (forming Pomeron exchange) are implied. 



^^abs ^Q i^YiQ inclusive structure function F2. Indeed, as we shall describe below, the application 
of the AGK cutting rules El E| gives^ 
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is the contribution to the diffractive structure function Fg (integrated 
over X]p) which originates from a perturbative component of the Pomeron of size l//i. Since the 
equality (jT} is valid for each component, yU, of the perturbative Pomeron, we can separate the 
screening corrections coming from low ^ < Qo (which are included in the input parameterisa- 
tions) from the absorptive effects at small distances (/x > Qq) which occur during the DGLAP 
evolution in the analysis of DIS data. Clearly, the inclusion of these absorptive effects will 
modify the parton distributions obtained from the DIS analysis. Not surprisingly, we find that 
by accounting for these 'Glauber-type' shadowing corrections we enhance the small-x input 
gluon distribution. 

A brief insight into the equality (Q) may be obtained from Fig. ^ which shows the leading 
absorptive correction (the two-Pomeron exchange contribution) to F2. Applying the AGK 
cutting rules jHI to this contribution, we obtain the relative contributions of -|-1, —4, and +2 
according to whether neither Pomeron, one Pomeron, or both Pomerons are cut. The first 
contribution is just F2 ■, which enters with the same magnitude, but the opposite sign, in the 
total /S.F2^^ . Hence the equality shown in (^Q). 

First, we summarise the fit to the DDIS data, which is based on a purely perturbative QCD 
framework jT]. The input Pomeron parton distributions were obtained from the lowest-order 
QCD diagrams for 7*p — > Xp \^. Recall that the crucial property was that the Pomeron wave 
function consists of a continuous number of components, i, each with its own size, 1/ ^i. The 
QCD DGLAP evolution of a component should start from its own scale /ij, provided that /i, 
is large enough for the perturbative evolution to be valid. Therefore, the expression for the 
diffractive structure function F2 contains an integral over the Pomeron size, or rather over 
the scale /i: 
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^Actually, to extract the leading-twist 'data' appropriate for a pure DGLAP fit of F2{xb,Q^) we liave to 
include the absorptive corrections AFI*^*^ integrated over fi'^ in the whole evolution interval from Qq to Q^. 



Here, xjp is the fraction of the proton's momentum transferred through the rapidity gap by 
the Pomeron, f3 = Xg/xp is the fraction of the Pomeron's momentum carried by the struck 
quark, Xg is the Bjorken x variable, and Q^ is the photon virtuahty. The subscript P on Fg p 
is to indicate that this is the perturbative contribution with /i > Qo — 1 GeV. The Pomeron 
structure function, F^(/3, Q^; /i^), is obtained by NLO DGLAP evolution up to Q^ of input 
quark singlet and gluon distributions of the Pomeron parameterised at a starting scale /i^. 
The (3 dependence of these input distributions is obtained by evaluating diagrams in which 
the virtual photon dissociates into either a quark-antiquark dipole or an effective gluon dipole, 
made up of a gluon and a compact qq pair.^ Bearing in mind the ii'-f actors arising from 
higher-order QCD corrections, the normalisations of the input distributions are treated as free 
parameters, Cg/c and Cg/c, to be determined in the fit to the DDIS data; see p. For the 
perturbative contribution, the Pomeron flux factor is given in terms of the gluon distribution 
of the proton,^ 
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where the subscript G indicates that the perturbative Pomeron is represented by two t-channel 
gluons in a colour singlet. 

In addition to the perturbative contribution to Fg given by (0), we also include a non- 
perturbative contribution (from scales /x < Qo), a twist-four contribution and a secondary 
Reggeon contribution, so that 

r^D{3) _ pD(3) pD(3) pZ)(3) pD(3) . ,x 

The detailed forms of these contributions are given in [Ij. In ;1] a good description of the recent 
preliminary ZEUS [3 EI and HI jlj DDIS data was obtained using a gluon distribution of the 
form Xjp g{x]p, ^'^) oc x]p^, with A ~ 0.17. However, the global parton analyses indicate that 
the gluon distribution, xg{x,fi'^), has a valence-like form at low scales /x ~ Qo ~ 1 GeV, which 
are dominant due to the l//i^ factor in the Pomeron flux factor Q, whereas the sea quark 
distribution, xS{x, /i^) = 2x[u{x, /i^) +d{x, iJ?) + s{x, /x^)], increases with decreasing x. For this 
reason, a sea quark-antiquark component of the Pomeron was introduced in ^, in addition to 
the usual two-gluon component. In analogy with Q, we therefore have a Pomeron flux factor, 
fjp=Si associated with this two-quark component, together with a GS interference contribution 

HI 

In this way, we are able to estimate the absorptive corrections, AF|'^''(xb, Q^; /i^), as a 
function of /i, from the perturbative component of F2 determined from the fit to the DDIS 
data. These absorptive terms can then be incorporated in the perturbative region of a new 
fit to the DIS data. Basically, the (negative) screening corrections have to be subtracted 
from the F2 data, before the DGLAP analysis is performed. At small x^, the F2 data, and 
consequentially the resulting parton distributions, are therefore appreciably enhanced. Notice 

^Earlier studies which first proposed the idea of determining the (3 dependence at the input scale, with 
subsequent DGLAP evolution, and which noted the consistency of this approach with DDIS data, were made 
in |S]. However, the DDIS data available then were less precise than at present. 

^To be precise, the flux factor should be written in terms of the skewed gluon distribution. At small xjp this 
gives rise to an overall constant factor ^0], i?^, which we absorb into the parameters Cq/Q and Cg/Q. 



that the original fit to the DDIS data jT] required knowledge of the gluon and sea quark 
distributions, xjp g{xjp, fi'^) and xjp S{xjp, ^'^), in the perturbative Pomeron fiux factors. Since 
the new DIS fit yields modified parton distributions, we therefore have to repeat the fit to 
the DDIS data. Fortunately, the successive iterations between the DDIS and DIS fits rapidly 
converge. 

We have described the fit to the DDIS data, so it remains to discuss the fit to the DIS 
data. Since we are primarily interested in the effect of absorptive corrections, it is sufficient 
to consider the description of the small Xg data. We therefore fit to the ZEUS jTTj and HI 
[12] F2{xs,Q'^) data with x^ < 0.01, 2 < Q^ < 500 GeV^ and W^ > 12.5 GeV^. We take 
MRST-like parametric forms JSI for the starting distributions at Ql = 1 GeV^: 

xg{x,Ql) = Ag x~^' {1 - xf-'^ {1 + eg^/x + -fgx) - A^ x~^- {1 - x)^^ , (5) 

xSix,Ql) = Asx-^'-il-xy-\l + esV^ + isx), (6) 

where the powers of the (1 — x) factors are taken from jl3,, together with the valence quark 
distributions, uy and dy, and A = d — u. The As, Aj, e^, A^, and 6- are taken as free 
parameters, with Ag determined from the momentum sum rule and 7^ initially fixed at zero. 
The charm and bottom quark contributions to F2 are evaluated in the light quark variable 
fiavour number scheme. The values of asiM"^) and the charm and bottom quark masses are 
taken to be the same as in the MRST2001 NLO parton set ^HI- Since we do not fit to DIS 
data with Xg > 0.01, we constrain the gluon and sea quark distributions to agree with the 
MRST2001 NLO parton set |T3j at x = 0.2. This is done by including the value of these parton 
distributions at x = 0.2 in the DIS fit with an error of 10%. To allow for the contribution of 
proton dissociation in ((T}, we take 

AFtixg, Q'; /i^) ^ -2F2^,,(x„ Q'; /i^), (7) 

where the factor 2 enhancement of the (elastic) proton contribution is estimated from the 
normalisation factors found in fitting to DDIS data [T|. The (elastic) proton contribution is 
obtained by normalising to the ZEUS LPS data |2], for which there is no proton dissociation. 
Integrating over /x^ from Ql up to Q^, we have 

F2U^B,Q'') = 6(0.1 -X,) f ' dxjp \f,'^J'\xjp,P,Q^) + F;^P{xjp,P,Q'')\ , (8) 

JxB 

where -Fg p is the leading- twist contribution and F^p is the twist-four contribution jT] . The 
upper cutoff of xp = 0.1 is necessary since the simple formula we have used for the absorptive 
corrections is invalid for large xjp (small rapidity gaps) where secondary Reggeon contributions 
become more important. The (negative) absorptive corrections, AFf^^, are subtracted from 
the F2 data, and the resulting leading-twist values fitted using NLO DGLAP evolution^ from 
the starting forms of (0) and ©. 

The dashed curves in Fig. El show the final parton distributions obtained after iterating 
between the fits to the DIS and DDIS data. The parameter values of this combined description 

*We use the QCDNUM program ^1] to perform the NLO DGLAP evolution and the minuit program JHl to 
find the optimal parameters. 





MRST-like input 


Xg = Xs = imposed 


F2 : xVd.o.f. 


1.09 


1.15 


^ 


10.1 (from mom. sum rule) 


0.82 (from mom. sum rule) 


A. 


-0.49 ±0.10 


(fixed) 


^9 


-1.2 ±0.1 


9.7 ±1.4 


l9 


(fixed) 


-15 ±2 


A_ 


(2.4 ±5.8) X 10-3 


(fixed) 


6_ 


0.74 ±0.30 


— 


As 


0.14 ±0.03 


0.56 ±0.04 


Xs 


0.30 ±0.02 


(fixed) 


es 


9.0 ±2.6 


4.0 ±1.2 


7s 


(fixed) 


-0.04 ±2.42 


F,""''^ : xVd.o.f. 


1.15 


1.29 


c,/G (GeV2) 


0.18 ±0.04 


0.37 ±0.03 


Cg/G (GeV2) 





(3.0 ±4.7) X 10-3 


CL/G (GeV^) 


0.074 ±0.032 


0.072 ±0.017 


c,/5 (GeV^) 


0.37 ±0.07 


0.032 ±0.007 


c,/5 (GeV2) 


1.14 ±0.07 


3.9 ±0.7 


Ci/5 (GeV2) 


0.027 ±0.033 





c,/NP (GeV-2) 


1.00 ±0.07 


1.38 ±0.05 


Cffi (GeV-2) 


6.5 ±0.5 


5.7 ±0.5 


iVz 


1.55 ±0.06 


1.53 ±0.06 


iV^ 


1.24 ±0.04 


1.20 ±0.04 



Table 1: The final parameter values of the two different simultaneous fits to the inclusive F2 
and F2 data measured by the ZEUS [21E1E] and HI |31^] Collaborations. The parameters 
for the F2 fit (279 data points) are defined in (0) and ©, while the parameters for the Fg fit 
(404 data points) are defined in pQ . Sample parton distributions from the two different analyses 
are shown in Figs. El and El 

of the DIS and DDIS data are given in the column headed 'MRST-like input' in Tabled The 
thin continuous curves in Fig. |21 show the parton distributions obtained from the fit before the 
absorptive corrections have been included; they are very similar to those from the MRST2001 
NLO parton set [121, with the input gluon distribution going negative at small x < 5 x lO^^. The 
dashed curves show that the inclusion of absorptive effects yield an input gluon distribution 
which is much less negative, whereas the input sea quark distribution is largely unaffected. 
Indeed, the absorptive effects crucially change the input gluon distribution for x < lO^^. They 
change the input sea quark distribution much less, due to the smaller colour charge of the quark 
and, phenomenologically, due to the fact that the quark distributions are measured directly by 
F2, whereas only scaling violations and NLO contributions constrain the gluon distribution. 
Thus small changes in the quark distributions can be accompanied by large changes in the 
gluon distribution. 

Note, from Tabled that the parameter A_ is consistent with zero. Indeed, repeating the fits 
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Figure 2: The gluon and sea quark distributions obtained from a NLO DGLAP fit to F2, before 
and after absorptive corrections have been included. The input at Qq = 1 GeV^ has been chosen 
to have 'MRST-hke' parametric forms, with an exphcit term included in the gluon distribution 
to allow it to go negative, see ^. 



with a fixed A_ = gives a description of the F2 data which is almost as good (x^/d.o.f. = 1.11). 
By contrast, without any absorptive corrections, the fit to F2 is much worse with a fixed A_ = 0, 
with a x^/d.o.f. = 1.57, compared to x^/d.o.f. = 1.15 if the input gluon distribution is allowed 
to go negative. 

Although the inclusion of absorptive corrections have enabled the DGLAP-based description 
to give a more physical small-x gluon distribution, they have not removed a long-standing 
puzzle of the behaviour of parton distributions at small x and low scales. That is, we still have 
a valence-like gluon distribution, whereas the sea quark distribution increases with decreasing 
X. That is, since the HERA F2 data have become available, we have had a 'Pomeron-like' sea 
quark distribution. Indeed, this feature has been present in all the parton analyses from GRV94 
[ini and MRS(A) [HI in 1994, up to the present CTEQ [IHl and MRST [19: global fits. On the 
other hand, according to Regge theory, the high energy (small x) behaviour of both gluons and 
sea quarks is controlled by the same right-most singularity in the complex angular momentum 
plane, and so we would expect 

A, = As, (9) 

where the Aj are defined in (0) and Q. If we impose such an equality on the Aj values, 
we obtain a very poor description of the F2 data. We have studied several possibilities of 
obtaining a satisfactory fit with this equality imposed, including saturation-motivated forms or 
including inverse transverse momentum ordering (which appears at NNLO), but none overcame 
the problem. The only modification which appears consistent with the data (and with the A^ = 
Xs equality) is the inclusion of power-like corrections. There may be higher-twist corrections 
due to the exchange of four gluons in colour antisymmetric states, which are not connected 
to F^ by the AGK cutting rules, and also more complicated higher-twist corrections caused 
by renormalons etc. Here we exploit the fact that such power-like corrections may slow down 
the DGLAP evolution at low Q^. Indeed, it has been argued [20] that such corrections must 
inhibit the growth of as and slow down the speed of evolution as Q^ decreases below about 
1 or 2 GeV^. At present, there is no precise formula to implement this effect. We therefore 
mimic it by shifting the scale ^J in F2{xb, Q"^) from Q^ to Q^ + m^, where m^ ^ 1 GeV^. To 
be consistent we must make the same shift in the Fg fit, so that, for example, becomes 

F^j'\xip,f3,Q^)= d^^ fip{xip;fi^ + m^)F^{f],Q^ + m^;fi^ + m^). (10) 

This simplified prescription enables us to obtain a satisfactory simultaneous description of the 
DIS and DDIS data, with the same asymptotic behaviour of the input gluon and sea quark 
distributions, as shown in Fig.|21 the corresponding parameter values are listed in the right-hand 
column of Table H 

A satisfactory simultaneous description of DDIS and DIS data has been obtained using the 
dipole saturation model |221I231- However, the description of the new more precise DDIS data 
using the BGK model ^21 is less good, with the model predictions tending to he slightly below 
the data, especially at low f3 [2j. Moreover, the DGLAP evolution of the Pomeron parton 
distributions is not accounted for. In the dipole approach, the best fit to DIS data also has a 
valence-like input gluon distribution [^ 1^ . This indicates that we need to account for the sea 
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Figure 3: The input gluon and sea quark distributions, (0) and ©, of our second analysis of 
DIS and DDIS data, in which the parameter A_ = and the equahty A^ = A5 (= 0) is imposed, 
as required by Regge theory. We now include 7j as free parameters in (0) and (0)- To obtain a 
satisfactory fit it is necessary to shift the scale in F2{xb, Q'^) from Q^ to Q'^ + 771"^, where m^ ~ 1 
GeV^. As before, the continuous and dashed curves show the parton distributions before and 
after absorptive corrections are included. We also show the 'shifted' parton distributions, at 
Q^ = {1 + 1) GeV\ 



quark contribution to the perturbative Pomeron flux factor in DDIS; indeed, tliis was one of 
the new ingredients of the analysis made in PP. Note that within dipole saturation models the 
sea quarks are generated solely from the gluon and therefore both have the same high-energy 
behaviour. In order to obtain a good fit to DIS data, the authors of |23 12^ were forced to shift 
the scale of the gluon distribution by /ig — 1 GeV^, the same value we used in (fTUj) . 

A phenomenological investigation of the effect of absorptive corrections from gluon recom- 
bination using the GLR-MQ J2S] equations was made in j^E] , before the advent of the HERA 
data. The input gluon and sea distributions were assumed to have a small-x behaviour of the 
form xg,xS ~ x^°'^ at an input scale of Qo = 4 GeV^. Since the small-x gluon distribution is 
now known to be valence-like at low Q^ from the HERA data, the shadowing corrections due 
to gluon recombination are correspondingly reduced, as found by MRST in 27j. At low Q^ the 
sea quarks are the dominant partons at small x, and hence sea quark recombination must be 
considered in addition to gluon recombination. 

Finally, a comment on why we consider partons at low scales. It might be argued that 
Q^ r^ 1 GeV^ is too low a scale to work in terms of quarks and gluons. (Recall that we only fit 
F2 data with Q^ > 2 GeV^.) However, we emphasise that Q^ ~ 1 GeV^ is the region where the 
description in terms of hadronic and quark-gluon degrees of freedom should be matched to each 
other. Therefore, we would like to obtain input parton distributions at Qo = 1 GeV^ which are 
consistent with Regge theory. An alternative approach is to adopt a hadronic description for 
Q^ ~ 1 GeV^ (see, for example, |2H|); however, this does not confront the issue. Note that within 
the operator product expansion (OPE), the leading-twist parton distributions are well-defined 
quantities even at low scales. Of course, at such low Q^, higher-order as corrections, power 
corrections and other non-perturbative effects are not negligible and need to be accounted for. 
Indeed, it was one of the goals of this paper to see if absorptive (and power) corrections could 
cure the anomalous behaviour of the gluon at low Q^ and small x. Note that the characteristic 
size of instantons, which are a typical example of the non-perturbative contribution, is about 
0.4 GeV^ (see, for example, [21]), and down to this scale it looks reasonable to work with quark 
and gluon degrees of freedom. The relevant hadronic (confinement) scale fih is smaller. It is 
driven by Aqcd and the constituent quark mass, that is, ;U^ ~ 0.1 GeV^. 

In summary, there is an outstanding dilemma in small-x DIS. Either, contrary to expecta- 
tions, the non-perturbative Pomeron does not couple to gluons, or DGLAP evolution is frozen 
at low Q^, perhaps by power corrections. Note, however, that in both scenarios we still have 
the puzzle that the secondary Reggeon couples more to gluons than to sea quarks. 
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